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1. Introduction 

2D string theory has thus far been most thoroughly studied in a specific background, 
defined by the standard gaussian model coupled to c=25 Liouville theory. An important 
open problem is to obtain an equally complete description of strings moving in other 2D 
backgrounds. 

One approach to this problem is based on perturbing the free action for two uncom- 
pactified real fields 0, X: 

^Liouville + ^Gaussian = J Z V§ W) 2 + £<W)) + J d 2 Z ^-{V X) 2 (1.1) 

by an operator e a *^C?£, where 0% are operators in the c=l gaussian model. In this paper 
we study the example of the Euclidean Sine-Gordon model coupled to 2D gravity: 



S = J d 2 zyf§ 
+ / d 2 z^/§ 



-^(VX) 2 + me w cos( P X/V2) 

Oil 



;i.2) 



Here g is some background metric. In flat space the Sine-Gordon model is not conformal 
for m 7^ 0. Coupling the Sine-Gordon model to gravity produces a nontrivial c = 26 
conformal field theory. General covariance (and hence conformal invariance) is maintained 
in the quantum theory for 7 = a/2, Q = a/8, £ = 7(1 — |p|/2). We generally use the 
notation, conventions, (and insights) of 0. The c = 1 model in 2D gravity is reviewed in 
§■ 

Correlation functions in the theory ( |1.2| ) will be defined by "conformal perturbation 
theory." That is, introducing the vertex operator 

V p = [ d 2 z^e^e ipX ^ (1.3) 



we define correlation functions at m 7^ by the series 



1 m ni+n2 
n li n 2>0 



The coefficients in the expansion (|1.4|) are calculated in the standard background with 
m = but \i 7^ 0. Recent results on the c=l matrix model have yielded a complete set of 
formulae for c=l correlators ||. In this paper we use these formulae to learn about the 
theory ( |1 . 2|) . Our main result is the phase diagram shown in figs. 2 and 3, and described 
in section four. Some physical interpretations of this diagram are proposed in section five. 
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2. Flows and phase transitions in 2D gravity 

In this section we review some (well-known) aspects of coupling constant flows in the 
c < 1 models coupled to 2D gravity. The discussion is meant to put the phase transitions 
discussed in sections four and five into perspective. 

The continuum approach to the c < 1 models begins with a u (p, q) theory" which is a 
tensor product of a Liouville theory with the minimal model M Pjq . Perturbations around 
this theory are defined by an action 

S = <Sp,q + T~r,kOr,k (2.1) 

r,k 

where the O r ^ are KPZ dressed operators in the Kac table, the latter being parametrized as 
in One of the couplings r r ^ must be nonzero to "set the scale," i.e., some operator must 
provide an infrared cutoff on the functional integral over surfaces. Correlation functions 
for nearby perturbed theories are defined by conformal perturbation expansions such as 
( |1.4|) . The nature of such expansions is difficult to analyze in the continuum theory since 
the coefficients are difficult to compute. Enter the matrix model. 

The solution of the continuum limit of the "(/-matrix model" indicates the existence 
of an infinite dimensional space of coupling constants, namely the space of real tuples 
{t rj /c}, 1 < r < q, < k, with all but finitely many t's = Q. There is a certain amount 
of evidence || that the continuum matrix model defined by the tuple i^fl = Sk q +r, P is 
identical to the theory S p ^ q . Just as the path integral for ( pTID is ill-defined if r r ^ = 0, 
the (p, q) string equation is singular for t Vj k — i^fl = 0. The mapping from t r ^ to T Vj k is 
nontrivial and has only been partially worked out in some special cases §. 

The advantage of the matrix model formulation is that the string equation and KP 
flow give a complete mathematical description of the crossover phenomena for coupling 
constant flow between the neighborhoods of two (p, q) "fixed points." Choosing some 
coupling, say x = t r0j k , to set the scale, the solution of the string equation u(x;t r ^) will 
be an analytic power series in the other couplings t Tt k- After proper identification of the 
Vs with the r's this power series should correspond to conformal perturbation theory. At 
genus zero the string equations reduce to algebraic equations for u. Hence, the power series 
in any coupling will in general have a finite radius of convergence. Therefore, if we make 
any coupling sufficiently large the physical description of the system must change. Let us 
consider two examples of this. 



2 



Example 1. The flow from Ising to pure gravity in the absence of a magnetic field is 
described by the equation 

u s + tu 2 = x (2.2) 

Letting x set the scale, for small t we may identify x with the cosmological constant and t 
with the coupling of the thermal operator e. The physical solution of (|2.2| ) is determined 
by the branch u = x 1 ^ 3 for x — > oo. For t < t c = (27X/4) 1 / 3 , u(x,t) is a convergent 
power series in t. In this regime the continuum description of the theory uses the action 
S = S'Liouviiie + Rising + t J e^e. Thus, for t < t c we are perturbing the Ising fixed point 
by a relevant operator. Beyond the radius of convergence the solution to ( j2.2|) must be 
expanded as a convergent power series in (x/t 3 ) 1 / 2 . We are now in the neighborhood of 
the (2, 3) fixed point and in the continuum theory we should describe this power series 
as a perturbation of the pure c = 26 Liouville theory by a certain irrelevant operator. 
Thus the action is S = SLiouvilie + t~ 3 / 2 J O where O is an operator in the c = 26 Liouville 
theory. Although the solution to ( |2.2| ) is analytic in the neighborhood of t = t c the a-model 
description of the physics changes. 

Example 2. It is easy to find examples of flows in coupling constants where there 
must be a true phase transition. Consider the string equation of the (2, 2m — 1) theories: 
u m + Yui>i ^i u% = x wri ere we let to = x set the scale. Consider the graph of the function 
f(u) = u m + Ysi>i ti 11,1 as a function of u, and denote the value of / at the local minimum 
with the largest value of u by h(U). If, as we change the couplings ti, h(U) crosses through 
x from below there will be a phase transition. In this case, if we simply analytically 
continue the specific heat around the branch point in complex t-space u will take complex 
values. 

These examples suggest a general idea, which is borne out by the results of this 
paper. There is a strong analogy between the coordinates on "theory-space," and 
weighted projective coordinates of complex manifold theory. First, gravitationally dressed 
operators depend on the Liouville zero mode only through a single exponential factor, hence 
the overall normalization of the r's can be changed by a shift of (p. Second, as we have 
remarked, the theory is singular if all the r r ^ = 0, reminiscent of the fact that the origin 
of <C n+1 is in no sense a point of <CP n . Third, in projective space regions in which a given 
coordinate can be scaled to one provide coordinate patches for the manifold. In r-space, 

1 Actually, this is only true for small t, and the results of || show that there is more to 
understand in this example. 
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different coordinate patches are defined by ietting different operators O r ^ set the scale. 
In example 1 above the change of expansion parameters from t/x 1 / 3 to x/t 3 is analogous 
to the change of coordinates between two coordinate patches of weighted projective space. 
In general, different coordinate patches correspond to different phases of the same theory. 

Let us apply some of these ideas to the Sine-Gordon model. Depending on the relative 
magnitudes of m and \x either the cosmological constant or the Sine-Gordon interaction 
will set the scale JIJ . If the cosmological constant sets the scale we expect that correlators 
will be expressed as power series in m 2 fi p ~ 2 . For example the genus zero partition function 
is Z = i / u 2 log / u + fi 2 f (m 2 fi p ~ 2 ) where f(z) has an analytic expansion around zero. On the 
other hand, if \i is small m sets the scale, and we expand in ^(m 2 ) _1 ^ 2_p - > . In section five 
below we will find that these expectations are in accord with matrix model calculations. 
We find a surprise in that there is a phase transition, possibly analogous to that of example 
2 above, for < p < 1 and for p > 2, while the model behaves much more like example 1 
above for 1 < p < 2. In section five we offer a physical picture that describes these phase 
transitions in terms of semiclassical field theory. 

Finally, we conclude with a few remarks on the relation between coupling-constant 
flow and renormalization group flow. A point where all but one r vanishes is rather like 
a fixed point of the renormalization group. At the fixed points we have a well-defined 
notion of matter central charge c x and bare matter-field dimensions. As in flat space, the 
operators perturbing away from the fixed points may be divided into relevant, marginal, 
and irrelevant. Since the Liouville field <p defines the local scale, the Liouville charge £ of 
the KPZ dressed operator e^&x is positive for relevant, zero for marginal, and negative 
for irrelevant operators. Thus, relevant operators grow in the infrared <fi — > +oo, etc. 
The classification of such operators is exactly the same as in flat space since by the KPZ 
formula and Seiberg's bound 0: f = § M - yjl + ^(A x - 1)^ gives f < for A x > 1 

and vice versa. In the Sine- Gordon theory for infinitesimal m the operator e ipX /^ 2 is 
relevant for p < 2, marginal for p = 2 and irrelevant for p > 2. More precisely, in flat space 
the renormalization group flow in the neighborhood of (m, p) = (0, 2) is given by the flow 
diagram of the Kosterlitz-Thouless model |7| . 

Ordinary renormalization group flow is dissipative ||. In gravity one may note, phe- 
nomenologically, that the flows by relevant perturbations of the g-matrix model between 
two fixed points always increases c x , while = c — 24A^ in always decreases. Con- 
versely, c^j-j- always increases under irrelevant perturbations. Thus, if we perturb by an 
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irrelevant operator and we find a phase transition we may expect that c^j has increased 
provided we have a phase transition to another surface theory. We will show in section four 
that one can perturb the Sine-Gordon model by an irrelevant operator p > 2 to obtain a 
phase transition. The extremely interesting question of whether this is a phase transition 
to a c^ff > 1 model remains open. 



3. The partition function for m^O. 

In this section we describe the result of a matrix model calculation of a one-point 



function of the form (|1.4| ). For technical reasons it is convenient to calculate the correlation 
functions of the vertex operators T± p = r( ._^ V± p , where p > here and hereafter. We 
define the coupling constant 

i r(-p) 
r(p) 



Using the calculational techniques of || and conformal perturbation theory we have found 
an explicit nonperturbative expression for the one-point function of the cosmological con- 
stant (7q) a = -^Z, where 

Z = ( e aT ^ +aT -^) (3.2) 

is the partition function. The somewhat complicated formula is given in equation (A. 6) of 
appendix A. In particular, defining a certain function, the "bounce factor" R q by: 

R q = ^-kly/| e W4 cos (| ( i +ifi _ M) ) r( i _ ifi+ (3 .3) 
we find that the amplitude 

A n (n, P ) = fi-^(%r^ p ) (3.4) 



may be expressed as a polynomial in the R q evaluated for g's at various integer multiples 
of p. See eq. (Al) for more detail. 
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3.1. String Perturbation Theory 

In order to study the partition function on a fixed topology we must expand the 
nonperturbative answer (A.l) in This may be obtained from the asymptotic expansion 
for the bounce factor: 



K p ~ exp 



^ (n + 1)! d/x 



1 + ^ + --- (3.5) 
A* 



where -00 denotes the expansion 

^ ~ E ^T/ 2 ^ 1 - 2 " 2fc+1 )4 - - - ^ (3-6) 

fc>i ^ 

$ is the digamma function and i?2fc are Bernoulli numbers. 

Substitution of (|3.5|) into equation (A.l) for A n (n,p) gives an asymptotic expansion 



of the form: 

An(»,p)~J2 u 2n-l + 2h A n(P) ( 3 ' 7 ) 
h>0 ^ 

The very statement of KPZ scaling, namely, that the above expansion begins at order 
2 *_i is somewhat miraculous from the point of view of the matrix model and implies the 
existence of nontrivial combinatorial identities on Bernoulli numbers. Nevertheless we may 
extract from ( |3.5|) and (A.l) the following facts about the correlation function: 

1. At each order of perturbation theory A^(p) is a polynomial in p of degree An— 2+Ah. 

2. p = is a zero of A^p) of order 2n. Indeed, as p — > 0: 

A n (^p)^p ln {^f n ~ 1 ReiP{\-in) 

(3.8) 



p 2n 



(2n-2)\ (2n)\ 



^Z 2 "" 1 24 / u 2n + 1 

3. The value p = 1 is a root of order n of 

4. Moreover, the value p = 2 is a root of order 1 for for h > 1, and in general for 
m a positive integer, p = m is a root of for h > ^(1 + n(m — 2)). 

Statement one is easily proved by examination of ( |3.5| ) . The expansion of the term in 
the exponent in powers p a /fx b has a maximum value of a — b for the term p 2 /[i. In || the 
amplitude (% Yli=i %>%) was shown to be a polynomial in bounce factors R p . It follows 
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that at genus h the amplitude is a polynomial in pi of degree 2k -2 + Ah. I A special case 
of this result is statement one. 

Statement two is easily proved from the low energy theorem in ||. As p — > we 
have T p — * pT = p-^- (except in the genus zero two-point function). Statement ( |3~8| ) 
immediately follows from the well-known value of the specific heat [pi . 

The proof of statements three and four is sketched in appendix B. 

3.2. The genus zero amplitude 

We now focus on the spherical topology. The above remarks show that A t ^ =0 (p) = 
(2n — 2)\p 2n (l—p) n Q n (p) where Q n is a polynomial of order n — 2 with Q n (0) = 1. Finding 
a formula for this polynomial has proved to be a rather difficult problem. Explicitly 
expanding the formula (A.l) we find, experimentally, the curious result 

n-2 

A£=°(p) =(2n - 2)\p 2n (l -p) n Y[(l-p/n) 

*=l (3.9) 

(T n T n v _ _ np-2n+2 , 2n ( -> _ xn T ( n ( 1 ~ P) + n ~ 2 ) 

'-W- ^ n ^ ^ p > r(n(l-p) + l) 

where r{ = 1 + i/n. We emphasize that this is a phenomenological formula, checked for 
1 < n < 11. We have made many unsuccessful attempts to prove ( |3.9|) for all n. i Two 
remarks might be useful to anyone else who tries: 

1. By explicit calculation, the special roots T{ are not roots of the genus one amplitude. @ 

2. The result can be summarized as a "Ward identity" similar to those which have been 
intensively studied recently in central New Jersey. Specifically, defining T p = pT p we may 
use the boundary-operator Ward-identity [1^] |13j to restate the result ( |01) as 

(t 7^ e t_V e } = {f?f\_ e f n e) (3-10) 

for < e < 1. 

In section four we will simply assume that ( |3.9|) holds for all n,p and explore the 
physical consequences. 



2 This confirms the observation of Q that at large energies the effective string coupling in the 

„2 



c = 1 model is g e f / ~ 2_ 

3 We thank R. Plesser for his participation in several of these efforts. 
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Thus the existence of these roots is reminiscent of the roots of the chromatic polynomial 



predicted by the Beraha conjecture [11]. 
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3.3. The special cases p = 1, 2 



Using discrete tachyon "Ward identities" recently derived in [fL3| we can give a much 
more complete description of the partition function and correlation functions for the special 
cases of a Sine-Gordon background with p = 1,2. 

At p = 1 the dependence on a is polynomial. Correlation functions at a ^ are easily 
related to correlation functions at a = 0. The general formula expressing this relation is 
somewhat long, so we simply quote a typical result 

(T T q T. q e^ +aT ^) = 2 F 1 (1 - q, 1 - q; 1; a 2 ^ l )iJ T q T. q ) (3.11) 

for q a positive integer. 

At p = 2 (when the matter perturbation is formally marginal) the T 2 Ward identity 
of |13| implies: 



Urn (t_ we f 2 "7? 2+e f[ % II t- qj ) = nl r( p f q " } (% {[ % \\ t_ qj ) (3.12) 

£ ^ i=l j = l i=l j = l 

where > 0, = Yllj = Q> an d Qi + Qj < 2 for all pairs It follows that 

correlation functions in this kinematic regime are given at a ^ by 

k i k i 

<n r ^ n T -^ eaT2+QT " 2 ) = (i-4a a r«<nr fc4 nr_ ft > (3.13) 

i=l j=l i=l i=l 

The correlation functions in other kinematic regimes will not be so simply related. For 
example, one can show 

(T q T_ q e aT2+aT -*) = (1 - 4a 2 )"« [l + Aa 2 q{q - 2)] (T q T_ q ) (3.14) 

for 2 < q < 4. 

Using the above methods one can directly derive the specific heat: 

(T r e aT2+aT - 2 ) = logfi - log(l - 4a 2 ) (3.15) 



which is, of course, in accord with (|3.9| ). An amusing, and perhaps important, feature of 
( P-15|) is that it is true to all orders of perturbation theory, since as noted in section (3.1) 
p = 2 is a root of the correlation functions at genus h > 1. 

We remark that the amplitudes ( |3.13| ) (|3.15| ) exhibit an interesting duality between 
the theory at a and at a = 1/ (4a) . This duality will be generalized below. 
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4. Phase Diagram in the a,p plane 

According to (|3J]) the genus zero specific heat is given by 



{T T e aT ?+ aT -i>} - (T T ) = -H(p; z) (4.1) 

where z = [i p ~ 2 a 2 p 2 (l — p) and we have defined an analytic function 

H(p;z) = Y JM 2 ^ , z» (4.2) 
^ n!r(n(l -p) + 1) V ! 

Some useful mathematical facts about the function H(p; z) are collected and proved in 
appendix C. In particular, H is a convergent power series in z for all real values of p 
with radius of convergence R c (p) given by (C.3) and plotted in fig. [I]. On the circle of 
convergence the series has one or two branch point singularities given by 

z c (p) = R c (p) p < 1 

zt (p) = -e ±l7ip R c (p) Kp<2 (4.3) 
z c (p) = -R c (p) p>2 

Moreover, the values of H(p; z) for \z\ > \z c \ may be related to values within the circle of 
convergence by connection formulae similar to those for the hypergeometric function. 

Using the above facts we may draw a phase diagram as shown in fig. [| There are six 
different regions: 

I. ) <p < 1, < /j p - 2 a 2 < R c (p)/(p 2 (l-p)) 

The specific heat is a power series in a 2 with real coefficients. There is a finite radius 
of convergence with singularity of the form (1 — z/zc) 1 ^ 2 . 

II. ) 1< p < 2, < n p ~ 2 a 2 < R c (p)/(p 2 (p - 1)) 

Again we have a real power series in a 2 , but the branch point has moved off into the 
complex plane, so there is no singularity as we approach the dotted line in fig. [| 

III. ) 2 < p < oo, < n p ~ 2 a 2 < R c (p)/(p 2 (p - 1)) 

There is again a finite radius of convergence with singularity (1 — z/zc) 1 / 2 since the 
branch point has moved back to the real axis. Of course, there is no singularity in passing 
between regions /, iT, iTi\ 

IV. ) <p < 1, /j p - 2 a 2 > R c (p)/(p 2 (l-p)) 
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The values in this region are defined by analytic continuation around the branch point. 
Explicitly, the connection formula (C.6) gives 

1 



{%%e aT - +aT ^) IV 



2-p 



+ 



1 



p-2 



log(tt 2 p 2 (l — p)) ± Z7T 

H 



p '._ e ±in/(p-2) 



(4.4) 



(aV(l-p)) 1 /(2-p)_ 

where p' = (2p — 3)/(p — 2). We have a power series in \x with complex coefficients. The ± 
sign depends on the sense in which we analytically continue around the branch point. 
V.) 1< p < 2, /j p - 2 a 2 > R c (p)/(p 2 (p - 1)) 

Now again applying the same connection formula, 



(T T e a ^ +aT -^ v 



' log(a 2 p 2 (p- !)) + ' 



H 



P 



A* 



(aVb-l)) 1/(2 " p) 



(4.5) 



2 — p P 
giving a power series in fx with real coefficients. 
VI.) 2 < p < oo, ^ p - 2 a 2 > R c (p)/(p 2 (p - 1)) 

In this region we use connection formula (C.8). The result is a power series with 
complex coefficients and expansion parameter 



(a 2 )- 1 /(p-D /U -(p-2)/(p-D 



(4.6) 



the expansion is analytic neither in a nor in p,. 

The case p = 2 requires special attention as discussed in section 3.3. 

Finally, we have used the coupling a which is natural from the matrix model. Changing 
variables a — > m using ( ^T|) the phase diagram in the p, m plane looks somewhat different, 
and is illustrated in fig. 0. 



5. Physical interpretation 

There are phase transitions when crossing the solid lines in fig. |2|. In this section we 
offer some qualitative physical interpretations of these transitions based on semiclassical 
analysis. The following considerations are only meant to be heuristic, and it would be 
interesting to make them more rigorous. 

The action ( |1.2j ) may be written as 

S =4^2 {/ + * + %m)] 

+ J d 2 zy^[^(VX) 2 +me^cos{pX)] 
10 



where for convenience we have rescaled and shifted 

m -> 47rm7 2 (2/^)^ 7 -»• 7 + log(/x/2) 

(5.2) 

For (5 = 2/7, £ = 1 we have a classical conformal field theory. Of course, quantum 
effects are strong at c = 1, but they can be summarized by the usual KPZ/DDK renormal- 
ization of parameters, so that 7 = y/2, Q = y/8, £ = 1 — p/2. Working semiclassically, the 
precise value of Q turns out to be unimportant so we will take the classical value Q = 2/7 
for simplicity. The physics depends sensitively on £ so we will leave this as a free parameter 
with < £ < 1. The equations of motion following from ( |5.1| ) are then 

R(e^g) + 1 + m£e" (1 -^ cospX = 
V 2 X + mpe^ sinpX = 

We restrict attention to the sphere with background metric: 

»=(ITH¥ * = 8 (B ' 4) 



(5.3) 



the constant solutions ("vacua") of (|5.3|) are given by (X n ,<p) where X n = ^(2n + 1), 



ti£ Z, and (f> solves 

8 + e*-fme^ = (5.5) 
This equation only has solutions for 

log(mO > (1 - 01og8 " ((1 " C)log(l - + ClogCJ (5-6) 

in which case the zero-mode potential looks like fig. (When there is no solution one 
must introduce Lagrange multipliers to fix the area of the surface. See |].) 

The existence of classical solutions for m larger than a critical value explains some 
features of the phase diagram of the previous section. In a phase where a solution exists 
we can expand around it and therefore we expect the partition function to be nonsingular 
for /u — ^ ]^] . This is in accord with the difference between regions /, iT and IV, V. 

When the condition ( |5.6|) is satisfied there are in fact two allowed constant curvatures, 
i.e., there are two solutions to (|5.5| ). For large values of m these are approximately 

e *.„(J*)i/* «1 

m i (5.7) 

~(m0 1/(W) > 1 
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The existence of different constant solutions, suggests the existence of solutions which in- 
terpolate between the different Sine-Gordon vacua X n and connect (fib to <fi s or (fib to (fib- 
We refer to such solutions as instantons, although this is an abuse of terminology. We 
may write a pair of ordinary nonlinear differential equations for radially symmetric solu- 
tions which, in the usual particle mechanics analogy ]n[] represent motion in the inverted 
potential —V((fi,X) as in fig. [5|. We have not proved that such solutions exist i but will 
assume that such solutions and their multi-instanton counterparts do exist, and proceed. 
The geometry associated to a multi-instanton configuration is of the form fig. [(]. 

We will now argue that the phase transition from IV — ► V may be thought of as a 
nucleation of a gas of instantons which unfreezes the classical vacuum (X n , (fib)- We may 
estimate the action of a single instanton for large m by the contribution of the potential 
energy: 

V((fi s ) - V(<fi b ) 



1 ^m^M-Dmi-o) (5 _< 



For on-shell configurations the kinetic energy should be comparable to ( |5.8|) . On the other 
hand, each instanton carries an entropy factor given by the relative areas of the large and 
small geometry 

~ 2(mO 1 /(«i-0) 8 - 1 /€ (5-9) 

Therefore, for £ < | and sufficiently large m the entropy will overwhelm the action and 
there will be a phase transition from nucleation of small bubbles, analogous to vortex 
unbinding in the Kosterlitz-Thouless transition. For £ > | this does not happen and the 
geometry remains frozen in a large sphere while the Sine-Gordon field remains frozen in 
its minimum. Our calculation is only meant to be a qualitative guide. The critical value 
£ c = | in fact coincides with the exact quantum result p = 1 of the previous section, but 
this is probably a coincidence. 

It is interesting to compare with the phase transition described in P| ||15|| . These 
papers show that an infinite chain of spacetime points sufficiently closely spaced act like a 
spacetime continuum in string theory. In our example, the chain of Sine-Gordon minima 
make X act like a discrete degree of freedom. The spacing between this chain of spacetime 
points is D = 2n(\/2/p) (we now revert to the conventions of the rest of the paper). 
Although we have not compactified X we may express D in terms of the minimal allowed 
radius of compactification. The transition at p = 1 corresponds to R = y/2, the self-dual 



The sphere metric complicates the usual analysis. 
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radius, while p = 2 is equivalent (under duality) to R = 2y/2, the Kosterlitz-Thouless 
radius . 

Let us finally consider region III. For p > 2, with the renormalized value of £ the 
zero-mode potential is unbounded from below in the ultraviolet. Nevertheless, the matrix 
model gives finite and real partition functions within a finite radius of convergence for m. 
(The matrix model definition of correlators of irrelevant operators implicitly makes a choice 
of regularization and subtraction of ultraviolet infinities.) The most mysterious transition 
is the the transition III — ► VI. Neither the Sine-Gordon coupling nor the cosmological 
constant can set the scale. This might be a signal of a first order phase transition to a 
"branched polymer" phase where the tachyon destroys the worldsheet 0. On the other 
hand, if it is a phase transition to a surface theory, according to the discussion of section 
two it could be a transition to a > 1 model. 



6. Nonperturbative Phase Transitions 



The effects of topology-change can drastically alter the nature of gravitational phase 
transitions. A good example of this is the difference in the nature of nonperturbative 
and perturbative phase transitions in the flow from the Yang-Lee edge singularity to pure 
gravity ||17|| . Analyzing the effects of topology-change in the phase transitions of Sine- 
Gordon gravity requires an understanding of the nature of the singularities of the full 
nonperturbative formula for the partition function given by equation (A.Q) below. 

Since the partition function is expressed in terms of a determinant: 



(T e aT ^ +aT -^) - (T ) = ilogDet 



(1 + £)(1 + E* 



(6.1) 



we must ask when the operator E, (defined in (A. 7)) ceases to be traceclass. From (A. 7) 
we may estimate 



< 



n>l 



(fi p a 



Rap I I R 



n>l 



a+6=n 

n(C(p)fi p a : 



bp I 



(a)! 2 (b\y 



(6.2) 



[niy 



where C(p) is a positive constant. For large n, R np behaves like (n!) p so for any < p < 2 
and small enough \a\, || E ||< 1 so the determinant exists and defines an analytic function 
of a 2 and p. 
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At p = 2 the sum in (0^2) behaves like a power series and remains convergent. More 
precisely, if we let p = 2 — is for s real we have for large n: 

\3lA < I e -^(*)+^i^- 2 ™(/( s ))- (6.3) 

Til 7T n 

where 

f(s) =4(l + s 2 /4)e^ |s| - s9(s) 

and 6(s) = tan" 1 (s/2). Thus if 4Ca 2 < 1 there is a finite interval — s c < s < s c along the 
Re(p) = 2 axis where the series is absolutely convergent. Thus there exists an analytic 
continuation of the determinant to the domain Re(p) > 2, where we are perturbing the 
Sine-Gordon theory by an irrelevant operator. 

We have just shown that, nonperturbatively, there are analogs of the regions /, II, III 
of section four. It is more difficult to see if the radius of convergence in [ck| will be finite. 
We expect that at fixed p there is a finite radius of convergence in | a \ . As we increase | a \ , 
1 + £ probably develops a left or right zero mode and the determinant has a logarithmic 
singularity, although we have not proven this. It is easy to show that for sufficiently large 
\a\, || E || > 1 so there is no reason for the determinant to be nonsingular. It would be 
interesting to understand the singularities better and to have a physical picture of how the 
phase transitions are modified by topology-change. 



7. Future Directions 

There are several projects which would extend the present work: 

1. Of course, it is important to prove ( |3.9| )! The recent results of [D5| are an important 
step in this direction. 

2. We also skipped over some hard analysis in section five, regarding the existence of 
instanton solutions. 

3. In section four we used the connection formulae (C.4), (C.6), (C.8). These relate different 
backgrounds via the action of the discrete group S3 and are thus reminiscent of target 
space duality and mirror symmetry. It would be very interesting to see if these symmetries 
survive in other correlation functions. 

4. We would like to have an equally complete understanding of the amplitudes at genus 
one which are implicitly contained in (A.l). These would be most useful for understanding 
better the nature of the transition JiT — > VI. If the transition is due to tachyon dominance 
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that should become apparent in the behavior of the genus one amplitudes. Unfortunately, 
we have not managed to recognize any special pattern in the first few amplitudes. 

5. The simple result (|3.15| ) deserves to be understood better. Naively, at p = 2 we have 
a tensor product of Liouville and Sine-Gordon theories, but for m 7^ 0, cos(v / 2A) is not 
exactly marginal so this is an illusion 0. Indeed ( |3.15| ) is not a product, but a sum of 
functions of (I and a 2 . 

6. It would be interesting to interpret the minimal models as restricted Sine- Gordon 
theories and relate the above results more directly to the c < 1 models. 

The present paper also touches on some deeper issues. The matrix model defines 
finite integrated correlators even for the irrelevant operators. How does it choose the finite 
parts? A natural guess for the underlying principle is the symmetry of the theory. 
This implies a corresponding symmetry of the continuum Liouville x matter system. 
Perhaps the finite parts are chosen based on the principle that the Ward identities 
must be maintained. 

A second issue is the spacetime interpretation of the Euclidean theory. If, for example, 
we compactify X when a = then we calculate the free energy of a string at temperature 
1/R. What happens when a / 0? The Euclidean Hamiltonian has now nontrivial .in- 
dependence. Should we interpret the calculations in terms of nonequilibrium statistical 
mechanics? 

One may also ask about the Minkowskian analog of the above results and the corre- 
sponding Minkowskian spacetime interpretation. Some of the relevant issues are discussed 
in 0. 

Finally, we may ask the evident question: Do analogous phase transitions exist in 
more realistic theories of gravity? 
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Appendix A. Derivation of the formula for the correlation functions 

We derive the formula for (To(T p T- p ) n ) by considering the e — > 0+ limit of the corre- 
lator (T e (T p _ e / n ) n (T- P ) n ) using the graphical rules of 0. A short calculation yields 

A n (^,p) = Li~ np (T T p n T™ p ) 

n , s k k k k d* (A 1 \ 

=.(-irM ! E^E(n^-n«» wnSi^ 

fc=l a,-,b ? - i=l i=l i=l 



where the sum runs over all partitions n = a± + b\ + ■ ■ ■ + a k + b k with ai, b\ > such that 
the denominator of 

C(a 1 ,b 1 ,a 2 ,b 2 ,...,a k ,b k ) = - — v — v , — — j-r -, — v — r (A.2) 

(ai + 0i)(6i + a 2 )(a 2 + t> 2 ) • • • (a fc + b k ){b k + ai) 

is nonzero. R p is the "bounce factor" of given by 

R p = ^~ P \/f el7r/4cOS (|^ + Z/U " P) ) r( ^ " ^ + P) (A - 3) 

for p > 0. 

The expression (|A1| ) can be written more succinctly by introducing an algebra of 1- 
dimensional projection operators V a ,b for a, b G Z + , not both zero, satisfying the relations 

/ai ' &l/a2 ' 62 (6i+a 2 )(6 2 + a!) /ai ' b2 1 J 

For such an algebra we may write the factor C as 

Such projection operators may be explicitly constructed as operators on a Hilbert space. 
Let \z a ), a = 0, 1, ... be an ON basis. Define \u>b) = J2 a ( a + ^) _1 | ;2 a) so that (zal^fc) = 
(a + 6)" 1 . Then V a ,b = (a + b)\z a ){w b \. 

Using the projection operators V a ,b the full partition function can be nicely expressed 
as a determinant of an operator S: 



(T e aT - +aT --) - (T ) = ilogDet 

where 



(1 + E)(1 + E 



(A.6) 



n>l a+b=n 
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Remarks: 

1. It might be an interesting exercise to obtain the above determinant directly from 
the fermion determinant in the original free-fermion formulation of the theory. 

2. Nonperturbatively the series vanishes identically at p = 0. This follows immediately 
since R P R* is real for p = and any real momenta p, q. 

3. In a similar way one can write slightly more complicated formulae for the three 
point function (ToT q T- q exp(aT p + aT- p )). 



Appendix B. Integer roots of the amplitudes 

The basic idea of the proof is very simple. To all orders of perturbation theory the 
bounce factor can be replaced by 

R _ u -P>p/2 r (l-^+P) 

- ^ (B.l) 

For p = n G Z + , R p becomes a polynomial in 1/p so that Qk(p = n) = for k > n. On 
the other hand, by KPZ scaling, we know the power of 1/n for the leading term in any 
amplitude. If this power exceeds the order of the relevant polynomials then we may prove 
vanishing theorems. For example, if we put p = m in ( |A.1[ ), then the expression must be 
of the form l/// nm times a polynomial in p. By KPZ scaling the genus h contribution goes 
like ~ 1/ n 2n - 1 + 2h so that the appropriate contribution must vanish for h > | + n(m — 2)/2. 

The proof that p = 1 is an n th order zero is much more tedious but uses the same idea. 
Having proved p = 1 is a root we take a derivative with respect to p of ( |A.1| ) . Plugging in 



p = 1 and using properties of gamma and polygamma functions we show that -§p-A n has 
an expansion in 1/p terminating at 1//U n . The proof then proceeds inductively and the 
inductive step fails when we consider (-^) n A n . 

We may note parenthetically that by the above reasoning any c = 1 amplitude with 
integral external momenta vanishes at sufficiently large orders of perturbation theory. This 
supports the general idea that special tachyons are associated with topological field theory. 
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Appendix C. Properties of the function H(p;z) 

In this appendix we prove some useful facts about the function 

_I>(2-p)) „ n 

n>l 



< n\T(n(l-p) + 1) 

1 ^ r(n(2-p))r(n(p-l)) . 
- 2^ n sin(n7rp)(-z) 



n>l 



71! 



(c.i; 



n> 



r(n(p-l)) 
' n!r(n(p - 2) + 1) 1 ' 



The ratio of gamma functions behaves at large n like 

n _3/2 (exp[(2 -p)log(2 - p) - (1 - p)log(l - p)])' 
n~ s/2 (-l) n sin(nnp)(exp[(2 - p)log(2 - p) + (p - l)log{p - 1)])' 



n 



-3/2 



{-l) n {exp[{p - 2)log(p - 2) - (p - - 1)])' 



< p < 1 
1< p < 2 

2 < p < oo 



(C.2) 



showing that the series defining H (p; z) converges absolutely for 



\z\ < R c (p) = exp 



(p — 2)log|p — 2| — (p — l)log|p — 1| 



(C.3) 



Note that from comparing the first and third lines we have the first connection formula: 



H(p; z) = -H(3-p;-z) 



(C.4) 



We can define analytic continuations of the function H (p; z) using various integral 
representations. 

Our first integral representation is the Mellin-Barnes representation 



1 

2ni 



r(s(2- P )) 
r( s (i-p) + i) 



r(- s )(-z) s ds 



(C.5) 



where z ^ 1R + and we use the standard branch of the logarithm. The integral over s 
converges absolutely for all such z if p < 1 and converges for \arg(z)\ > n(p — 1) for 
1 < p < 2. For \z\ < \z c (p)\ we can close the s integral in the right half- plane to obtain 
the series in the first line of (|C.1| ). If \z\ > \z c (p) \ then we can close in the left half-plane, 
thus proving the connection formula 



h(p-z) = -J_iog(-z) + -^-if(p';-(-^r 1/(2 - p) ; 



p 



2 -p 
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(C.6) 



for p' -2= -V 

Our second integral representation is the Mellin-Barnes integral 



-1 

2ni 



r( 8 (p-i)) 

r(a(p-2) + l 



-r(-s)z s ds 



(C.7) 



where z ^ Hi - . For p > 2 this converges absolutely for all such 2 and defines an analytic 
continuation of H. For 1 < p < 2 the integral converges absolutely for |ar<7(z)| < 7r(p — 1). 
By closing into the right half-plane we obtain the series in the third line of ( |C.1|) and by 
closing in the left half-plane we obtain the third connection formula: 



H(p;z) 



1 — p 



1 \ogz + -^H(p'-z- 1/ip - l) ) 



p — 1 



(Ci 



where p' — 1 = —^r. 
1 p—i 

Our third integral representation is derived from the middle series in ( |C1| ) using the 
integral representation of the Beta function. The result is 



Z ~Hz H(yP; ^ 



z sin irp 



n 



dt- 



-ty 



(l + ze in n 2 -p(l - ty- 1 ) (l + ze- i7T H 2 -P(l - t)^ 1 ) 

(C.9) 

This integral always converges at the endpoints t = 0, 1 for 1 < p < 2 and defines an 
analytic continuation in z for these values of p. The existence of singularities can be 
examined by looking for pinching of the contour. In this way it is easy to check that there 
are no singularities as z increases from zero to infinity through real values. 

In general H (p; z) does not seem to be expressible in terms of standard special func- 
tions, although at some special values we can write H more explicitly: 



H(0;z) = -lo£ 



1 + y/T=Iz~ 



dt 



2 

+ -7T 

3 



1 2 1 27t 2 
2-^1(0' o-> o! ~ T- 



3' 3 2 



dt 

T 



(C.IO) 



H{l-z) = -\og{l-z) 



H(-;z) =21og 



z z z 
1 + T + 2 



if(2;z)=log(l + z) 
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When p is rational, H (p; z) can be written in terms of generalized hypergeometric functions 
a Fi,. In general for fixed p the branch point singularity in z is a square root singularity, 
except at p = 1,2 where we have a logarithmic singularity. 

Finally we note that the set of transformations of the (p, z) plane defined in 
( P-4|) ( |C.6|) (|C.8|) defines an action of the permutation group S3 on this plane. 
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Figure Captions 



Fig. 1. Radius of convergence as a function of p. 

Fig. 2. The phase diagram in the a 2 vs p plane. The solid lines indicate lines across 
which there are phase transitions. 

Fig. 3. The phase diagram using the normalization standard for vertex operators. There 
are now infinitely many separated regions where we have a transition of the type 
III -> VI. 

Fig. 4. The zero mode potential V(0, X) = 8(j) + + me^ cospX in the case when there 
are constant postive curvature solutions to the equations of motion. 

Fig. 5. Particle motion of (X(r), 4>(r)) for a proposed field configuration connecting large 
and small geometries. 

Fig. 6. A multi-instanton configuration. Many small instantons of scale size r 2 ~ e^ s 
join onto a large sphere of scale size r 2 ~ e^ b . 
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